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For a bi-partite quantum system denned in a finite di- 
mensional Hilbert-space we investigate in what sense entan- 
glement change and interactions imply each other. For this 
purpose we introduce an entanglement-operator, which is 
then shown to represent a non-conserved property for any 
bi-partite system and any type of interaction. This gen- 
eral relation does not exclude the existence of special ini- 
tial product states, for which the entanglement remains small 
over some period of time, despite interactions. For this case 
we derive an approximation to the full Schrodinger-equation, 
which allows the treatment of the composite systems in terms 
of product states. The induced error is estimated. In this 
factorization-approximation one subsystem appears as an ef- 
fective potential for the other. A pertinent example is the 
Jaynes-Cummings model, which then reduces to the semi- 
classical rotating wave approximation. 



I. INTRODUCTION 

During the last decades entanglement has been investi- 
gated under various aspects. The famous EPR-paradox, 
for example, has led to a discussion of the most basic 
principles of quantum physics The Gedanken ex- 
periment based on "Schrodinger's Cat" may be seen as 
an attempt to challenge the consistency of quantum me- 
chanics: It has been argued that this situation could only 
be understood by allowing for entanglement between the 
atom and the cat which, on the other hand, should be 
considered a classical object 0|. But, by definition, a 
classical object cannot become entangled with any other 
system. 

Since it has been shown, that quantum algorithms have 
the potential to outperform corresponding classical com- 
puting H H |Q, considerable efforts have been made 
to implement gates like the so-called quantum controlled 
NOT-gate (QCNOT). Performing a QCNOT generically 
results in preparing an entangled state. Meanwhile var- 
ious experimental schemes to prepare entangled states 
have been developed [|| , @ • 

In all these approaches entanglement has been in the very 
center of interest. The question was always cither how 
to interprete the state of two systems being entangled, or 
how to deliberately produce entanglement and detect it, 
once it has been produced. 

Rather neglected seems to have been the question of en- 



tanglement as an unavoidable "waste product" of quan- 
tum mechanical dynamics. Little attention has been paid 
to the fact that it cannot be taken for granted that any 
two interacting systems will remain in a product state, 
even if they have been in one in the begining jll|] 
This means that there is always the possibility for them 
to entangle. And if they are entangled, it is impossible 
to assign two separate wavefunctions to the subsystems. 
Nevertheless this is typically done in standard " textbook 
level" quantum mechanics: The particle in a box, e.g., 
is always described by a wavefunction although it defi- 
nitely interacts with the box that neccesarily consists of a 
many particle-quantum-system itself and therefore could 
become entangled with it. There is no discussion of the 
electron going through the double slit being possibly en- 
tangled with the material defining the slit itself. 
But since these approximations typically lead to excel- 
lent results, it should be possible to point out why. In 
which situations is it reasonable to neglect entanglement 
and treat whole complicated systems as effective poten- 
tials for another quantum-system? 

Apart from the rather academic desire to undestand the 
basis of this " classical limit" , there is also a good practi- 
cal reason to adress such questions. 
An important prerequisit of all quantum computer de- 
signs suggested so far is the possibility of so called lo- 
cal unitary transformations. These should be performed 
selectively on each effective spin (q-bit) through poten- 
tials that are supposed to be controllable in time jl5). 
But again, in reality, those potentials can only be imple- 
mented by means of other complicated quantum-systems 
that could possibly entangle with those spins: this would 
inevitabcly lead to decoherence. But quantum comput- 
ers need to be coherent. In that sense the problem of 
entanglement through interaction (as required by exter- 
nal control) could even challenge the implementation of 
any real quantum computer. 

Our paper is organized as follows: We first specify a theo- 
rem relating some purity measure P (as an entanglement 
test) to inter-subsystem interactions (Sect. [nj). For the 
proof of this theorem (Sect. |n|) we proceed as follows: 
Starting from the Von-Neumann-equation, which de- 
scribes the dynamics of the density operator, we pro- 
ceed by inserting an expansion of the density operator 
into this equation. The result is an equation only in 
terms of the expansion coefficients that has exactly the 
form of the Schrodinger-Equation and will therefore be 
called " quasi-Schrodinger-equation" . It is now possible 
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to define a linear operator in the space of those coeffi- 
cients which has an expectation value equal to P, and 
will therefore be called " purity operator" . Since the dy- 
namics of those coefficients are controlled by an equation 
that is formally identical with the Schrodinger-equation 
(including a "quasi-Hamiltonian"), it is possible to re- 
duce the question of P being conserved or not, to the 
problem whether the commutator of the purity operator 
and quasi-Hamiltonian will vanish or not. 
Thus, the mathematical scheme used here is essentially 
the same as used in standard quantum mechanics to iden- 
tify conserved quantities. Only the space of the state 
vector and the interpretation of the considered quanti- 
ties, are different. 

The last step will be to show, that the above commutator 
becomes nonzero whenever the full Hamiltonian involves 
any kind of interaction. 

However, even in the presence of interactions the sys- 
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tern may remain "almost" unentangled. In Sect 
we use our quasi-Schrodinger formulation to derive the 
factorization-approximation with its effective potentials 
for this case. In Sect. [v| the induced error is estimated 
to lowest order. In Sect. [v| we apply the results to the 
Jaynes-Cummings-model. 



II. THEOREMS 

There is a still ongoing debate on entanglement mea- 
sures ||. A lot of propositions have been made, but it 
seems still rather difficult to introduce a general entan- 
glement measure that satisfies all conditions that have 
been imposed on such a measure and, at the same time, 
is applicable for any number of subsystems and any case 
(pure and mixed states of the whole system). And it 
seems even more difficult to construct a measure in such 
a way that it could actually be calculated (or measured!) 
for reasonably complicated situations. 
Fortunately, it is possible to introduce a simple measure 
under specific conditions: If the state of the whole system 
is a pure state, and the full system is being regarded as 
divided into two subsystems, a convenient entanglement 
measure is 1 — P, where 



P = Tr / {(/) 2 }=Tr // {(p 
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(2.1) 



(here p 1 , p 11 are the reduced density matrices of the cor- 
responding subsystems) 

Entanglement between to subsystems originating from 
unitary quantum evolution, can only result from inter- 
actions []. If two systems do not interact they can be 
treated without even taking the other one into account. 



1 Non- unitary transformations can do without direct interac- 
tions: This phenomenon has become known as entanglement 
swapping JL6| . 



So, if they are both in pure states at the begining, which 
means they are in a product state with respect to the 
whole system, they will remain so forever under these 
conditions. 

One may ask now whether two systems that interact 
might remain entanglement-free, depending for example 
on the kind of systems that interact, or on the kind of 
interaction that is considered. 

Concerning this question we are aware of only rather 
vague statements in the literature. A typical formula- 
tion due to d'Espagnat reads: 
Theorem A: 

"In general it is impossible to describe systems that in- 
teracted in the past by separate wavefunctions" O]. 
But does this always have to be the case? To address 
this problem we will prove the following theorem for fi- 
nite discrete Hilbert-spaces: 
Theorem B: 

"There exists no interaction what so ever between ar- 
bitrary systems, such that the entanglement measure 
(1 — P) remains conserved". 

This theorem does not imply that there cannot be initial 
states, starting from which the sytem might remain in a 
product state, though it can be shown that those states, 
if they exist at all, only play a negligible role in typical 
larger systems. But it definitely means, that there must 
be initial states that lead to entanglement, even between 
the particle and the box-system, or between the electron 
and the slit-system. 

Further consequences of theorem B can most conve- 
niently be assessed from an approximation-scheme that 
is valid as long as the systems remain approximately un- 
entangled, as will be shown in Sect. |lV|, ^ and Sect. VI. 



III. PROOF OF THEOREM B 

A. Basis operators 

The operators into which the density matrix will be 
expanded here, are products of the generators of the re- 
spective SU(n) groups, where one set of generators cor- 
responds to one subsystem || . 

The basis operators for the different subsystems are de- 
fined in the following way: 



-v 



1 \ v 



i = 
i ^ 



(3.1) 



where, v denotes the index of the subsystem, n v the num- 
ber of levels of subsystem v (dimension) and AJ" the gen- 
erator of the SU(n u ) group, i = 1, 2 . . . n 2 — 1. 
The basis operators of the full system are defined as 
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dyadic products of the basis operators of the subsystems, 

N 



(3.2) 



where N is the number of subsystems. 
Thus any basis operator is defined by a sequence of N in- 
dices i v (abbreviated as i), each index specifying, which 
generator i should be applied to the corresponding sub- 
system v. The operators constructed according to these 
rules form a complete and orthonormal set in the follow- 
ing sense: 

N 

(3.3) 

where A is an arbitrary operator. 

Representing the density matrix in the case of only two 
subsystems as 



(3.4) 



where the index i corresponds to subsystem I and the 
index j to subsystem II, we find 



'mi 



E qio Q* p = nu E q 



(3.5) 



The objects we are going to examine aer thus specified 
in terms of their expansion coefficients. 



H*^ = Tr {H [Qj 7 Q, Tl ] }* = Tr j (h [Qj, Qa])^ (3.9) 

We can even define a formal bracket notation: For this 
purpose we re-arrange the multiple indices j as a sim- 
ple index s and introduce a set of real orthogonal basis 
vectors 



with 



such that 



s s 



H ss > = (s\n\s') q s 
Inserting ( |3 . 1 2 ) into (3.8) yields: 



(s\q) 



ih-\q)=H\q) 
with the infinitesimal solution 



\q(dt)) = ( 1 + -Hdt ) \q(0)) 



(3.10) 
(3.11) 

(3.12) 
(3.13) 
(3.14) 



On the basis of this formal equivalence it is now possible 
to investigate the conservation of some quantity A in the 
space of the vectors \q) {(q\A\q) = ^ s s , q s A ss >q s >), by 
evaluating the commutator of A with TL. 



B. Quasi-Schrodinger-Equation 

The Von-Neumann-equation reads: 



Inserting the expansion (3.4) yields: 
d 



(3.6) 



(3.7) 



After multiplying by Q^, taking the trace and applying 
some trace theorems, we get: 

4™ = E% n Aj :=Tr{H[Qj,Q m ]} (3.8) 



This equation (hereafter called quasi-Schrodinger- 
equation) has evidently the Schrodinger-form. The hcr- 
miticity of TL is easily shown by examining the corre- 
sponding matrix elements, 



C. Purity Operator 

We first note that our purity measure, P, can be given 
the mathematical form of an expectation value: 



(q\V\q)=P 



(3.15) 



To find V , we go back to the explicit multi- index notation 
(for two subsystems), 



P ^ y q&P ss' q& ! — ^ y q%jPijuvqu 



(3.16) 



l,J,U,V 



If we require (see (|3.5[)) 

(q\V\q)^ niI J2qfo (3-17) 



it follows that 



(3.18) 



i. e., up to a normalization factor V is a projector, pro- 
jecting out the components of \q) that refer locally to 
subsystem I, those components that would read in 
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multi-index notation for a bi-partite system. Using the 
bracket notation V reads: 



(3.19) 



Adding the complementary projector that projects out 
all the other components, that is all those that do not 
refer locally to subsystem I, we can write the unity oper- 
ator as: 



(3.20) 



Using this representation Tt, which controls the complete 
dynamics, can be split up into: 



H = c 1 + n + w 



(3.21) 



where 



C 1 :=^£LlO(m 7 | 

l.rn 



ho 

1,3 

One easily convinces oneself, that 

[P,£ / ]=0 and [P,n]=0. (3.22) 
Finally, the commutator [P, W] reads: 

[V,W] = n 2 n J2 (W™jW){3 R \ - Wj^Kw 1 ]) . 

v>,3 

(3.23) 

We now convince ourselves that not all matrix-elements 
yV w k can be identically zero in the presence of interac- 
tions. For this purpose we define "interaction" opera- 
tionally: Two subsystems are said to interact, if the dy- 
namics of one subsystem, at least for some initial state 
(but generically for any initial state) depend on the state 
of the other subsystem. 

The state of subsystem I at time dt is completely dis- 
cribed by the projection 



\q\dt)) :=£|«/)<«/|g(d0). 



(3.24) 



Substituting \q(dt)) according to the evolution equation 
(3.14) with Ti given by (3.21) we obtain 



\q\dt)) = |/(0)) (3.25) 
+ Ut\Y J LiJw I ){m I \ q (0))+Y J W w] \w I ){ ] %(Q)) 



As only the components (j R \q) carry information about 
subsystem II not all 1/F-matrix-elements can be zero for 



interacting subsystems. 

Since the two terms that are summed over in ( 3.23 ) ob- 
viously belong to different (off diagonal) parts of Tt, they 
cannot cancel each other. Thus we conclude, 



(3.26) 



for any sort of interaction, which completes the proof of 
theorem B. 



IV. FACTORIZATION-APPROXIMATION AND 
EFFECTIVE POTENTIALS 

As we have shown, there is, a priori, no reason to as- 
sume that two interacting systems remain in pure states. 
Nevertheless, as we will argue now, it is a reasonable ap- 
proximation to treat one subsystem as if it was in a pure 
state, and the other one as an effective potential for the 
former, as long as 1 — P remains small over the period of 
time under investigation. 



A. Factorization-approximation 

We start with the quasi-Schrodinger-equation ( |3.§| ) for 
two subsystems, 



ih- Qij 



EM 



H 



qir, 



(4.1) 



and consider only those equations which refer to the two 
subsystems locally, not to the correlations. We thus need 
to study: 



I : ihy/nJJ—q l0 
dt 



■ 5> \h 



II : ih^/nj-^qoj = - ^ Tr lit 

l,m \ 



vr 

Al q 



qim (4.2) 



I in 



qir 



i.e. we have to examine the relation between "local" 
{ctiOiQoj) and "global" (q^) coefficients. We define a ten- 
sor M in the following way: 



Mij :={Qi< 



(Qi)(Q 3 - 



or, according to ( |3.4|) and (3.3) 



Mij = qtj - s/rwiijqioqoj 



(4.3) 



(4.4) 



Obviously, for a product state all Mij 's have to vanish. If 
the two subsystems get correlated (entangled) , the My 's 
will take on nonzero values, so that 



(4.5) 
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can be considered as an alternative entanglement mea- 
sure. It is indeed possible to derive a relation between (3 
and P M: 



[3<l-P 2 



(4.6) 



I f w e now solve (4.4) for <fc and insert the result into 



(4.2) we get for subsystem I: 



ih^/nJJ—q i0 
at 

-£Tr{tf 

l.m \ 



y/2' 



(4.7) 

{y/njnjjq m q 0m + M im ) 



Using (p|) and (lO) we find 



nin u 



\ ~> 2 2 



1lQ%r, 



p 2 j2 M ^ <1 - p2 



(4.8) 



so that it seems reasonable to neglect the M; m -term and 
only keep the ^/ninnqio gom-term in (4J), as long as P 
stays close to 1. The error that occurs if this approxi- 
mation is used over a period of time instead of the true 
Schrodinger-equation, will be estimated later. 
With the above approximation we get after performing 
the sum over m: 



ihy/nJJ—q i0 
at 



H 



A t Inn ~ 



qia 



(4.9) 



The latter is a quasi-Schrodinger-equation for subsystem 
I with a quasi-Hamiltonian depending on the momentary 
local state of subsystem II. It is easy to check that the 
quasi-Hamiltonian of this equation, although depending 
on time, will always be hermitian. This means that in 
this approximation ^/nJiJ2i Qio — P IS a conserved quan- 
tity. Both subsystems will remain in pure states and un- 
entangled if the initial state was a product state. 
Performing the sum over I and taking the partial trace 
with respect to subsystem II yields: 



ih^/nji—qio = 
at 



-Trj I Tin {tip 11 } 



v/2' 



(4.10) 



Since pi and f>2 always represent pure states in this case, 
we can use separate wavefunctions to re-express lokal ex- 
pectation values, such as: 



Tr 



ii 



{hp 11 } = ^ n \m n ). (4.ii) 

Attention should be paid to the fact that this expectation 
value is still an operator with respect to subsystem I, if 
H contains interactions. 

Following the standard trace theorems, we get: 



(4.12) 



Multiplying the i-th equation by -^t, summing over i, 



exploiting the completeness of the ^|'s and using (3.5) 
leads to: 



lh dt Pl 



(4.13) 



This equation is of the same form as the Von-Neumann- 
equation. Since the reduced density operators pi, pn rep- 
resent pure states we can change to the corresponding 
Schrodinger-equation, without loss of generality: 

iftj^ J > = (^I^I^W) (4-14) 

After splitting the Hamiltonian H into local (Li, Ljj) 
and interaction (W) parts and adding the corresponding 
equation of subsystem II, we find the complete system of 
equations: 



Ui 



+ {i> n \W\i> 11 )) \^) (4.15) 



(ip u \L n \i> 
B. Gauge 



Since the overall phases of both subsystems are ar- 
bitrary and can be chosen freely, a method, which is 
completely analogous to the Lorentz gauge in classical 
electrodynamics, can be used to simplify these equations 
even further. 

Applying the substitutions 



and choosing the phases as 



ja u (t)\ill\ ._ 



ii \ 



{^LiW) 



f- Mi(t) = 
hau(t) = 



(4.16) 



(4.17) 



transforms the system of equations into: 



= (l u 



h tf\w\tf)\ W 1 ) 



(4.18) 



This is a coupled set of nonlinear first-order differential 
equations which, although it does not necessarily cre- 
ate unitary dymamics, keeps the absolute values of the 
wavefunctions of both subsystems fixed. Instead of the 
(n/ • nn) dimensions of the exact treatment, this ap- 
proximation has only (rii + nn) dimensions (number of 
equations). 

Each of these equations can be considered as an "ordi- 
nary" Schrodinger-equation of one system, in which the 
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influence of the other one appears as an effective poten- 
tial. This approach may be termed " quasi-classical" with 
respect to interactions. It underlies, e.g. the potential for 
the particle in a box. 

Similar equations are well-known in the theory of many 
particle systems (for example the Hartree-equation [[l2| ) . 
But in those cases they are basically justified by their 
success. We are now going to derive a criterion for their 
applicability. 



V. ERROR ESTIMATION 

We will now examine the deviation of the solution in 
factorization-approximation |0(i)) from the solution of 
the full Schrodinger-equation \ip(t)), under the condition 
that |0(O)) = \ip{0)) ■ We will expand this deviation in 
terms of an effective interaction which will allow us to 
connect the resulting error with the purity P. 
The Schrodinger-equation reads: 



%n^\il)) = [L, +L,, + \Y 



(5.1) 



With \(j>) :— {(j) 1 ) (g) |</> // ) the factorization-approximation 
reads: 

ih^\<j>) = (Lj + L u {tf\W\<?) + tf^WW 1 ) + a) |0) 

(5.2) 

where a denotes an insignificant overall phase, and 

|0) := |0 7 ) ® |0") (5.3) 

Introducing the deviation \9) as 

\l/>) =: \6) + |0), (5.4) 

and the abbreviation 

V :=W - (0 / |M / |0 / ) - (<t>"\W\<t> n ) - a (5.5) 

we can write a differential equation for this deviation \9) 
just in terms of |0) and \0) as 



d 



dt ih 
Defining 

\ff) :=e^ 
we can rewrite as 



H\6) + V\4>) 



(5.6) 



(5.7) 



The exponential may be written as 

e^=J2\MmMt)\=Y,\MV))({<t>n(t)\ + {0n(t)\) 
n n 

(5.9) 



Here the |0„(O)) = |?/>„(0)) are chosen to form a complete 
orthonormal set. The |0 n (t)) as solutions in factorization 
approximation form a complete orthonormal set at all 
times if they did so at t — 0, so they will be use d as a 
basis. Since they all obey (5.8) we get, using (5.9): 



n 

(5.10) 

Transforming back to the unprimed quantities, we find 
\6 t (t)) = (5.11) 

4 E + Km /' ^(t')\ + (e n (t')\) v\Mt'))dt> 



By iterating this equation we can produce an expansion 
of \0i(t)} in terms of time integrals over matrix elments 
of the effective interaction V. As long as those matrix el- 
ments remain small, a truncation scheme can be applied, 
which gives to first order: 



!**(*)> = 4 E [ (<t>n(t')\VMt'))dt' (5. 

We return now to the double-index notation 



12) 



\<ki) := |0() ® |0f ), (<f>M') =&i>i, WW) = Sj>i 

(5.13) 

and take |V'oo)(|0oo)) as the solution(approximation) un- 
der consideration. We use the |0y)'s as a basis to write 
the full solution 



IV'oo) = /l - E l%l 2 l^oo) +5Z%|0ij), #oo := 

(5.14) 

Its overlap of this solution with the solution in 
factorization-approximation ("fidelity") is exactly given 
by 



i-Ei 



(5.15) 



hi 



On the other hand, the squareroot of the purity of either 
subsystem of the state IV'oo) to first order in |%| 2 is given 
by 



Vp = i- l%l 



(5.16) 



(5.8) 9. 



Obviously the two quantities are the same except for the 
Vs and the Ooj's. Computing the Oij's from equation 
( ^12| ) yields: 



\l (<0y|W|0oo)- (00^1^1000)^0- (5.17) 



(0io|W^|0oo)^oj + aSioSoj) dt' 
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Choosing a := (^ool^l^oo) to satisfy the definition 

#00 : = from ( 5.14 ) we get 



so that 



0io(*)=O %(*)=0 



|<V>oo|<M| 2 = V^ 



(5.18) 
(5.19) 



as long as the small deviation of \f~P from 1 is domi- 
nated by the lowest order term in the effective inte ract ion 
strength, as it is the case in the example in Sect. VI. 
This means that no "local" errors are generated in first 
order, any deviation occuring due to the factorization- 
approximation leads to entanglement and therefore re- 
duces purity. 



VI. APPLICATION TO THE 
JAYNES-CUMMINGS-MODEL 

A. The model 

The factorization-approximation will now be applied 
to the Jaynes-Cummings-model. This choice is moti- 
vated by the fact that this model can be solved exactly, 
hence it can be tested whether the approximation is re- 
ally applicable. Furthermore, both methods of treating 
this system, the fully quantum mechanical and the semi- 
classical one, are well known, thus the result of apply- 
ing the factorization-approximation can easily be inter- 
preted. 

The Jaynes-Cummings-model Jl3| describes a spin in 
a magnetic field, interacting with some monochromatic 
electromagnetic field. A good example is a typical NMR 
experiment of any sort. 

The model in rotating wave approximation is defined by 
the Hamiltonian: 

H = Bg/j,-a z + Hu)(a^a + — ) + hj(a^a^ + aa + ), (6.1) 

where B denotes the magnetic field, g the gyromagnetic 
relation, /i Bohr's magneton, a z the operator of the z- 
component of the spin, 7 the coupling constant, a + ,a~ 
are the creation and anihilation operators of the spin sys- 
tem, u> is the freqency of the electromgnetic field and 
a) , a are the creation and anihilation operators of the 
electromgnetic field. 

The second term describes the monochromatic electro- 
magnetic field, the first term the spin in the magnetic 
field and the third term their mutual interaction. 



B. Application of the factorization-approximation 



(6.2) 



(6.3) 



= [Bg^* z + h 7 ((<f> L \a)\ct> L )a- + (cj ) L \a\^ L )a+) ) \<, s ) 

= Uuj Ua+^j +ft 7 «0 s |a-|^)at + (0 s |a+|0 s )a)^ \4> L ). 

where S indicates the spin system and L the system of the 
electromagnetic field. The absolute values of the expec- 
tation values {4> s \a~ (j) S ),{4> s \a + 4> s ) are always limited 
by: 



\(d> a \&-<^)\, \{<P S \a + 4> S )\<l 



(6.4) 



The influence of the coupling terms on the evolution of 
the electromagnetic field is thus negligible, if 7 <C u> 
and/or the system of the electromagnetic field is in a 
highly exited state. Thus, for example, a coherent state 
I a) with a large parameter a is a valid solution to the 
equation controling the dynamics of the electromagnetic 
field. 

Inserting this solution into the equation for the spin sys- 
tems leads to: 



(6.5) 



= \Bg^a z + fry ((a\a + \a)a- + (a\a\a)cr + )^ \<f) S ) 

= yBgp,^a z + fvyV2\a\ (cos(ujt)(j x + sm(ujt)a y )^ \4> s ), 

This is exactly the semi-classical approach, which of- 
ten produces exellent results, such as the correct Rabi- 
frequency etc. 

It remains to be shown, that it was indeed justified to 
use the factorization-approximation for some time r. 



C. Estimation of P 

If the spin and electromagnetic field are in resonance, 
an exact solution of the Jaynes-Cummings-model can be 
found: 



1 



\n+) : =-^(|0>®|n) + |l)®|n-l) 



(6.6) 



V2 



(|0) ® |n) - |1) ® |n - 1) 



Applying the scheme (4.18|) to this system yields: 



which satisfies the following eigenvalue equations: 

H\n+) = (hun + h-yy/^)\n + ) 
H\n~) — (hum — hjy/n)\n~). 



(6.7) 
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Considering the initial product state 



In the same limit we also find: 



®|o), 



(6.8) 



we then find for the components of the Bloch vector of 
the spin system: 

(a x ® 1) = -Im (KAn+ie- lut ) (sin (7 (V^ + v^+T) t) 



v / (« + i)~2Vn >/»- + ~°- ( 6 - 15 ) 

For a large enough a ( |6.9| ) therefore simplifies to: 

(6.16) 



sin(a;i) |A„| 2 sin (27\/ni) 

ri 

cos(wi) I A„ 1 2 sin (27V™*) 



- sin (7 (Vn - \/n + 1) t)) 
1) = 5^ Re (A* n A n+1 e~ luJt ) (sin (7 (Vn + V^TT) t) 



E 



-4, 



; (2 7v ^t) . 



in (7 (Vn- Vn + 1) t)) 



sin 

E 



I A, 



; (2y/nt) 



In this approximation the factors depending on u>t de- 
scribe the Larmor-precession, while the sum terms give 
rise to the Rabi-precession. 
(6-9fexpanding the argument of the Rabi-precession arround 
the mean value of the Poisson-distribution yields: 



From these components we find P as: 
Tr s {p%} =P=±(l + (a x ) 2 + (a y y 



(6.10) 



For the initial state we again choose a coherent state (for 
simplicity with a real parameter a), which means: 



2 a" 



(6.11) 



Now we want to simplify the terms in (6.9). Consider: 

v«+l 



A n A n 



+1 



-\a\< 




(6.12) 



The first factor describes a Poisson-distribution which is 
characterized by the mean value a 2 and the standard de- 
viation a. Since the main weight of such a distribution is 
concentrated near its mean value (98% within 3 standard 
deviations) , it is a reasonable approximation to keep only 
contributions from this range. The second factor can be 
estimated using its respective value at the boundaries of 
this range: 



< 



a 



< 



TL re i 



a ■ 



1) 

(6.13) 



where n re i denotes the n's from within this relevant 
range. If we now again only consider highly exited co- 
herent states, that means the limit a — * 00, the upper 
as well as the lower bound converges against 1. Hence in 
this limit: 



yj(n re i + 1) 



1. 



(6.14) 



(6.17) 



7* [2a + — (n - a ) 
\ a 



2a 3 



(n-a 2 ) 2 + O d (n-a 2 ) 



Since we are only considering contributions from the rel- 
evant region, that is contributions with \n — a 2 \ < a 
in the limit a — > 00, we can neglect any term beyond 
second-order and get: 



(n — a 2 ) 



(6.18) 



Using this approximation and so me ad dition theorems 
we get on the right hand side of ( 6.10 ) fast oscil lating 
terms with a frequency given by the first term in ( 6.18j ) 
and slow oscillating terms with a frequency given by the 
second term in (6.18). Since the slow oscillating terms 
from the relevant region oscillate with a frequency of the 
order of 7 regardless of the actual a, we expand those in 
time to second order, which will be a good approximation 
as long as jt remains small: 



P = 1 



(6.19) 



This means that in the limit of high intensities of the elec- 
tromagnetic field, the "coherence time" is independent of 
the intensity itself and only depends on the strength of 
the coupling. If that is weak enough, the factorization 
approximation will be valid for a considerable time. 



VII. SUMMARY AND CONCLUSIONS 

We have shown that bi-partite systems with mutual 
interactions, that remain unentangled forever, no matter 
which initial product state was chosen, cannot exist. 
One might argue that this is a rather weak statement: 
In large systems with many degrees of freedom there 
might be a large number of initial states, starting from 
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which the subsystems would remain entanglemcntfrce, 
even though there must also be some initial states that 
lead to entanglement. 

But the same could be argued, for example, for a classi- 
cal Hamilton-system, which has no radial symmetry, say: 
Although angular momentum is not conserved in such a 
situation in general, there might be special trajectories 
for which it is conserved. But one probably would not 
claim that these trajectories are of importance consid- 
ering the space of all possible trajectories, especially in 
systems with many degrees of freedom. 
The same holds true for the entanglement conservation: 
there might be special initial states that prevent the sys- 
tems from getting entangled, but it can be shown that 
their relative weight decreases with the system size. 
Thus it is, strictly speaking, unjustified to describe a par- 
ticle in a box, which is part of an interacting quantum 
system, by a wave-function. 

It is, nevertheless, a potentially very good approximation 
to describe interacting systems by pure individual wave- 
functions, if the entanglement remains small during the 
time of observation. The underlying equation can then 
be based on the factorization approximation; this equa- 
tion has been derived, together with a criterion for its 
validity. 

For a spin interacting with a monochromatic electro- 
magnetic field, we have shown that the criterion is ful- 
filled for high enough field intensities and for a time that 
only depends inversely on the coupling strength. Apply- 
ing the factorization approximation then transforms the 
full quantum mechanical Jaynes-Cummings-model into a 
semi-classical one, in which the field acts as an external 
potential. 
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